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Abstract We investigate the effect of exotic matter in 
particular, hyperon matter on neutron star properties 
such as equation of state (EoS), mass-radius relation- 
ship and bulk viscosity. Here we construct equations of 
state within the framework of a relativistic field theo- 
retical model. As hyperons are produced abundantly in 
dense matter, hyperon-hyperon interaction becomes im- 
portant and is included in this model. Hyperon-hyperon 
interaction gives rise to a softer EoS which results in a 
smaller maximum mass neutron star compared with the 
case without the interaction. Next we compute the coef- 
ficient of bulk viscosity and the corresponding damping 
time scale due to the non-leptonic weak process including 
A hyperons. Further, we investigate the role of the bulk 
viscosity on gravitational radiation driven r-mode insta- 
bility in a neutron star of given mass and temperature 
and find that the instability is effectively suppressed. 

Keywords Neutron stars, dense matter, r-mode 
oscillation 
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1 Introduction 

The investigation of spin frequencies from burst oscilla- 
tions of eleven nuclear-powered millisecond pulsars showed 
that the spin distribution h a d an upper limit a t 730 
Hz (jChakrabartv et all 120031 IChakrabartvl 120051 ). The 
fastest rotating neutro n star discovered re cently has a 
spin frequency 716 Hz ( Hessels et al.l l2006h . In this re- 
spect, the study of r-modes in rotating neutron stars has 
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generated great interest in understanding the absence of 
very fast rotating neutron stars in nature. The r-modes 
are subject to Chandrasekhar-Friedman-Schutz gravita- 
tional radiation instability i n rapidly rotating neutron 
stars (lAnderssonlll998ll2003HFriedman fe Morsinklll998t 



Lindblom, Owen & Morsink 1998t lAndersson. Kokkotas fc Schutd 



Il999t IStereioulasI 12003) which may play an important 
role in regulating spins of young neutron stars as well as 
old, accreting neutron stars in low mass x-ray binaries 
(LMXBs). Another aspect of the r-mode instability is 

that t his could be a possible source for gravitational radi- 

ation (lAnderssonl2003|:lGondek-Rosinska. Gourgoulhon fc Haensell 
l2003tlNavvar fc OwerfeOQq ) which may shed light on the 
interior of a neutron star. 



It was realised that the r-mode instability could be ef- 
fectively suppressed by bulk viscosity due to exotic mat- 
ter in neutron star interior when the compact star cools 
down to a temperature ~ 10 9 K. The coefficient of bulk 
viscosity due to non-leptonic weak processes involving 
exotic matter such as hype ron and quark matter was cal- 
culated by several authors ([J ones l200lHLindblom fc Ower] 



2002 



2004; 



Haensel. Levenfish fc Yakovlevl2002tlvan Dalen fc Dieperinkl 
Draeo. Lavagno fc Pagfiard2005tlMadsenlll992l . l2000[) . 



It was noted that relaxation times of the non-leptonic 
processes were within a few orders magnitude of the pe- 
riod of perturbations and gave rise to large bulk viscosity 
coefficient. This led to the complete suppression of the 
r-mode instability. 



In this paper, we investigate the effect of hyperon 
matter including hyperon-hyperon interaction on bulk 
viscosity coefficient and the r-mode stability. In Sec. 2 of 
this paper, we discuss the model used to calculate equa- 
tions of state, bulk viscosity coefficient and the corre- 
sponding time scale. Results of our calculation are anal- 
ysed in Sec. 3. Section 4 gives the outlook. 
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2 Model 

2.1 Equation of State 

We describe the /3-equilibrated and charge neutral hy- 
peron matter within the framework of a relativistic field 
theoretical model where baryon-baryon interaction is me- 
diated by the exchange of scalar and vector mesons and 
hyperon-hyperon interaction is taken into account by two 
strange mesons - scalar f n (hereafter denoted as a* ) and 
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Kf, 
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vector 4> (|Schaffner et al lll993t lMi shustin fe Schaffnerll99 ' 
The Lagrangian density for hyperon matter including 
non-interacting electrons and muons is written as 

C-B = ^s(«7m^ M ~ m B + g<jB<? + 9a'B<J* 
B 

- 9uiBl^ - g^Bl^ ~ g P Bl^B ■ P^B 

+ \ (d^a - mla 2 ) - U(a) 



where kp B , Jb, Hb and ns are Fermi momentum, spin 
and isospin projection and number density of baryon B 
respectively. The last terms in energy density and pres- 
sure are due to leptons. The effective baryon mass is 
defined as m* B — tiib — g a B& — 9a* b&* ■ The charge neu- 
trality and chemical equilibrium conditions are given by 
Q — J2b Ibtib - n e - n M = and //, = bifi n - q^ where 
^t n , [i e and [ii are respectively the chemical potentials of 
neutrons, electrons and baryon i and bi and qi are baryon 
and electric charge of ba ryon i. The chemical potential 
of ba ryons B is given by (jChatterjee fc Bandyopadhyavl 
2006) [i B = (k 2 F B +m B 2 ) 1 l 2 +g ujB u +g < t > B4 l o+hBgpBPm- 
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The Dirac bispinor Wb represents isospin multiplets for 
baryons B—N, A, E and E. Here vacuum baryon mass 
is tub, and isospin operator ts while w M „, §pv, and p M „ 
are field strength tensors. Neutrons and protons do not 
coupl e with a* and <j> mesons. T he scalar self-interaction 
term (jBoguta fc Bodmer1ll977t ) is given by 



U(cr) = ^g 2 i 3 + \g^ 



(2) 



We perform the calculation in the mean field approxima- 
tion. The total energy density and pressure are respec- 
tively given by 
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and 
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In this calculatio n we adopt nucleon-meson co upling 
constants of GM set (jGlendenning fc Moszkowsklll99lh . 
The coupling constants are obtained by reproducing prop- 
erties of saturated nuclear matter such as binding en- 
ergy — 16.3 MeV, saturation density n = 0.153 /m~ 3 , 
asymmetry energy 32.5 MeV, effective nucleon mass 0.78 
and incompressibility 240 MeV. The hadronic masses 
used in this calculation are m a = 550 MeV, m u = 783 
MeV, Trip = 770 MeV and m N = 938 MeV. On the 
other hand, hyperon- vector meson coupling constants are 
determined using SU(6) symmetry of the quark model 



(IMishustin fc Schaffnerll99ftlDover fc7^alll984l : [Schaffner et all 
Il994h and the scalar a meson coupling to hyperons is cal- 
culated from hyperon potential dept hs in normal nuclear 



matter such as U^(n ) = -30 MeV dDover fc Gal|[l98 



Chrien fc Dover! 1989h Uff(nn) = -18 MeV dFukuda et al.1 



19981 : iKhaustov et all [2000) and a repul sive potential 



depth for E hyperons (n„) = +3 MeV (iFriedman Gal fc Battvl 
11994 iBattv. Friedman fc Ga]|ll997l : iBart et al.lll999f ) as 
obtained from hypernuclei data. The hyperon-cr* cou- 
pling const ants are determined from double A hyper- 
nucle i data ( Schaffner et al ] l!993t IMishustin fc Schaffnen 
Il996h . 

We obtain the EoS solving equations of motion along 

with the expression for effective baryon ma ss and charge 

neutrality and beta equilibrium constraint s (iMishusti n fc Schaffner! 
119961 : IChatterjee fc Bandvopadhvavll2006h . 



2.2 Bulk viscosity coefficient, damping time scales and 
critical angular velocity 

Now we discuss the bulk viscosity coefficient in young 
neutron stars which cool down to temperatures ~ 10 9 — 
10 10 K after their births. As the system goes out of 
chemical equilibrium due to pressure and density vari- 
ations associated with the r-mode oscillations, micro- 
scopic reaction processes in particular non-leptonic weak 
interaction processes including exotic partic l es migh t re- 



particles 

store the equilibrium (iNavvar fc Owenll2006t [j oncs 2001 
Lindblom fc Owenl l2002t IChatteriee fc Bandvopadhva - 



2006b . Here we calculate the real part of bulk viscosity 
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Similarly, we compute other components of atij . It is to 
be noted here that nuclcons do not couple with strange 
mesons, i.e. g a * n = 0. 

The hyperon bulk viscosity damping timescale (th) 
contributes to the imaginary part of the r-mode fre- 
quency and is calculated in the following way dNavvar fc Owen 



1999) 



20061: iLindblom. Owen fc Morsinkll9 98: Lind blom. Mendell fc Owed 
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Th 



1 dE 
2E~dt 
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where E is the energy of the perturbation in the co- 
rotating frame of the fluid 



E : 



e(r)r dr 



(12) 



Fig. 1 Mass-radius relationship of non-rotating neutron 
stars including nucleons-only matter (blue solid line) and hy- 
peron matter with (dashed line) and without (magenta solid 
line) hyperon-hyperon interaction. 



The derivative of E with respect to time is given by 

dE '* 
dt 



-47T 



C(r)(|V -5w\ 2 ydr, 



(13) 



where we need to know the energy density e(r) and bulk 
viscosity £(r) profiles of a neutron star. Similarly, we es- 
coefficient (0 in terms of relaxation times of microscopic t imate the damping time scale (tu) corresponding to the 
processes (jhandau fc Lifshitdll999l : lLindblom fc Owenl2002t)i iillr viscosity due to the modified Urea process including 

P(loo - 7o)t 



1 + (U)T) 

where the difference of infinite (7no) and zero (70 ) fre- 



3ity due to trie modified Urea process : 
on ly nucleons usin g the bulk viscosity coefficient as given 
bv lSawverl ((Hi). It was noted that gravitational radi- 



quency adiabatic indices is g iven by (jLindblom fc Owenl 
l2002t (Navvar fc OweJxM) 

P dn n drib 

In the co-rotating frame, the angular velocity (u>) of 
(I, rn) r-mode is related to angular velocity (J?) of a ro- 
tating neutron star as ui — mppp, ^ ( Anderssonll2003l ) . In 

this calculation, for I — m — 2 r-mode, it is 10 
The relaxation time (r) for the non-lcptonic process 

n + p ^p + A , (7) 

is given bv ILindblom fc Owenl (|2002h 

- - 3 k F A (\M A \ )- 
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ation drives the r-modes unstable. So the gravitational 
radiation time scale (tgr) has a negative contribution 
to the imaginary part of the r-mode frequency. Further, 
the overall r-mode time scale (r r ) is defined as 

1111 

- = + — + — . (14) 

T r T GR T B T V 

The r-mode is stable below the critical angular velocity 

which is obtained from the solution of — = 0. The critical 

Tv- 
angular velocity depends both on the mass of the neutron 

star and its temperature. 



3 Results 



192vr 3 



n b 5x n 



(8) 



where kF A is the Fermi momentum for A hyperons, (\Ma \ 
is the angle averaged squared matrix element and 

S " (9) 



UbSXn 

with a 



01 An - OL n A + OLAA 



We obtain expressions for ctij from 



the baryon che mical potential and equations of mo tion 

for meson fields (|Chatterjee fc Bandyopadhyay|200i) . For neutrality condition. Due to the appearance of additional 



We p erform this calculation using the GM set (jGlendenning fc Moszkowsli 
I1991D as it has been described in sec 2. It is found that A 
hyperons appear first at 2.6n followed by E~ hyperons 
at 2.9no. However no £ hyperons appear in the system 
beca use of strongly repulsive i7-nuclear m atter interac- 
tion ( Chatteriee fc Bandvopadhvavll2006T ). Further, we 
find that after the appearance of negatively charged S 
hyperons, number densities of electrons and muons de- 
crease in the system. This is attributed to the charge 



example, 

d^A 



ttAr, 



9u]A9ujN 



dn n 



da 



da* 



Ik 



dn n 



(10) 



degrees of freedom in the form of hyperons, the EoS 
of hyperon matter becomes softer compared with that 
of nucleons-only matter. On the other hand, the inter- 
play between the attractive a* field and repulsive <j> field 
which becomes dominant with increasing density, makes 
the EoS softer initially and stiffer at higher densities than 
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Fig. 2 Bulk viscosity coefficient as a function of radial dis- 
tance for a non-rotating star (solid line) and as a function of 
equatorial distance for a rotating star (dashed line) of mass 
1.6 Mq including hyperon matter with hyperon-hyperon in- 
teraction. 



Fig. 3 Critical angular velocity as a function temperature 
for a rotating star of mass I.6M0 including hyperon matter 
with hyperon-hyperon interaction. 



the EoS without hyperon-hyperon interaction. This be- 
haviour of the EoS is reflected in the calculation of max- 
imum mass and the corresponding radius of the neu- 
tron star. In Fig. 1, we plot the mass-radius relation- 
ship of non-rotating neutron stars calculated by solv- 
ing Tolman-Oppenheimer-Volkoff equation for equations 
of state including nucleons-only matter (blue solid line) 
and hyperon matter with (dashed line) and without (ma- 
genta solid line) hyperon-hyperon interaction. Here the 
highest maximum mass corresponds to the stiffest EoS 
of nucleons-only matter. The maximum masses corre- 
sponding to the EoS with and without hyperon-hyperon 
interaction are 1.64 M Q and 1.69 M Q respectively. Using 
the rotating neutron star (RN S) model of Stergioulas 
( Stergioulas fc Friedmanlll995t ). we find that the maxi- 
mum masses corresponding to the EoS with and without 
hyperon-hyperon interaction are respectively 1.95 Mq 
and 2.00 M . 

Next we discuss the hyperon bulk viscosity due to the 
non-leptonic process given by Eq. (7). The bulk viscosity 
coefficient is dependent on the EoS as it is evident from 
the expression in Eq. (5). The temperature dependence 
of the bulk viscosity enters through the relaxation time 
as given by Eq. (8) . The bulk viscosity coefficient reaches 
a maximum an d then drops with increasing baryon d en- 
sity as found bv lChatteriee fc Bandvopadhvavl (|2006l ). It 
has been also noted that the bulk viscosity increases as 
temperature decreases. Therefore, the large value of ( 
might be effective in suppressing r-mode instability as 
neutron stars cool down to a few times 10 9 K. 

The neutron star which was rotating rapidly at its 
birth, slows down due to gravitational wave emission. 
A rotating neutron star has less central energy density 
than its non-rotating counterpart because of the centrifu- 
gal force. Consequently, hyperon thresholds are sensitive 



to rotation periods of compact stars. In the calculation 
of critical angular velocity, we consider a rotating neu- 
tron star of mass 1.6 Mq and the corresponding central 
baryon density is 3.9no which is well above the thresh- 
old of A hyperons. This star is rotating at an angular 
velocity 2952 s -1 . It slowed down from its Keplerian an- 
gular velocity 5600 s _1 when the central baryon density 
was below the A hyperon threshold. The calculation of 
critical angular velocity requires the knowledge of en- 
ergy density and bulk viscosity profiles in compact stars 
as it is evident from Eqs. (11)-(14). In Fig. 2, the hy- 
peron bulk viscosity profile (dashed line) of the rotat- 
ing neutron star of mass 1.6 Mq is shown as a function 
of equatorial distance along with that of a non-rotating 
(solid line) star of same mass. We note that the peak of 
the profile shifts towards the center in the case of the 
rotating neutron star. Earlier calculation of critical an- 
gular velocity was performed using energy density and 
bulk vi scosity profiles of non-r otating neutron stars. Re- 
cently, [Na^iya7&^Owen| ( 20061 ) studied the effect of ro- 
tation on critical angular velocity. We calculate the crit- 
ical angular velocity using energy density and bulk vis- 
cosity profiles of the rotating star. In Fig. 3 the critical 
angular velocity is plotted versus temperature. We find 
that the r-mode instability is suppressed by hyperon bulk 
viscosity. Here we do not consider the impact of super- 
fluidity of b aryons on bulk viscosity. However, several 
calculations (iHaensel. Levenfish fc Yakovlevll2000L l200ll 
l2002tlNavvar fc Owenll2006D showed that superfluidity of 
baryons might suppress the bulk viscosity and the damp- 
ing of the r-mode instability. 
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4 Outlook 

Other forms of matter such as K~ condensed matter 
and quark matter may appear and compete with the A 
hyperon threshold in neutron star interior. In an earlier 
calculation we showed that the condensate of K~ mesons 
appeared at ~ twice normal nuclear matter density fol- 
lowed by the appearance of A hyperons (jBanik fc Bandvopai 
1200 ll) . It would be worth investigating how the bulk vis- 
cosity coefficient due to the process (7) is modified in this 
situation. Also, we are studying whether the non-leptonic 
process involving K~ condensate n ^ p + K~ could give 
rise to bulk viscosity coefficient as large as that of the 
non-lepton ic process (7) and damp the r-mode instability 
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